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<n: 

The Kronig-Penney model with random Dirac potentials on the lattice Z has critical 
energies at which the Lyapunov exponent vanishes and the density of states has a van Hove 
singularity. This leads to a non-trivial quantum diffusion even though the spectrum is 
■ known to be pure-point. 

a- 

1 Main results 



>• ; The Kronig-Penney model describes the motion of a particle in a one-dimensional system with 
singular potentials. We consider the case of a random 5-potential on the points of the lattice Z. 
p\j ■ The Hamiltionian is given by 

where u = (v n ) ne z are i.i.d. random variables. We assume that the distribution of the i> n 's has 
compact support and is nondegenerate, i. e. is not concentrated in a single point. The precise 
mathematical meaning of the ^-potential in ([I]) is recalled in Section [2] below. If the f n 's all 
have the same sign, the almost sure spectrum of H u consists of an infinite number of bands, 
while the spectrum contains all positive reals if is in the support of the distribution of the f n 's 
[KM] (see Section [2] for details). Under rather general assumptions it is known that the spectrum 
the one-dimensional random Schrodinger operators is pure point and that the eigenfunctions are 
exponentially localized. We are not aware of any publication in which this is actually proved 
for model ([I]), but the general methods [Stoj can be applied for this model (details will be given 
elsewhere, see, however, |HKK] where localization at the bottom of the bottom of the spectrum 
is proved). 

The main focus here is on the transport properties of the random Kronig-Penney Hamiltonian. 
We analyze the growth of the time-averaged q tb moment of the position operator X on L 2 (M): 

M q (T) = - dt e~% e %Hut \X\ q e- tHut q > . (2) 
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This operator is an integral operator and its integral kernel is denoted by (x\M q (T)\y). The 
diagonal entries (a\M q (T)\a) can be interpreted as the time-averaged moments of a wave packet 
initially localized in a Dirac state at a G R (which, of course, is not an element of Hilbert space). 
The following result shows that for sufficiently large q these moments grow with time T, which 
is certainly not the typical behavior inside the localization regime. 

Theorem 1 Suppose that E(u) 7^ and let a G R \ Z. For every a > there is a positive 
constant C a such that for T > 1 

(a\M q (T)\a) > C,^" 8 , q> h -. (3) 

The heuristics leading to a non-trivial lower bound ([3]) on the wavepacket spreading is as 
follows. We will show that the localization length (given by the inverse of the Lyapunov exponent) 
diverges at certain critical energies E\. At the same time, the density of states diverges at these 
energies so that there are many such states (see Theorem [2] below). The quantum motion in these 
states is more or less ballistic until it reaches the localization length. This decreasing fraction of 
delocalized states allows to prove the lower bound. Both these heuristics and the implementation 
of the idea in a proof are quite similar to the random polymer model |JSS] . There is, however, 
a crucial difference: the transfer matrices at the critical energies E\ are all Jordan blocks in 
the random Kronig-Penney model, while they are random rotations at the critical energies of a 
random polymer model. This leads to a completely different behavior of the Lyapunov exponent 
and ultimately also for the moments in ([3]). Indeed, in the random polymer model the behavior 
is (a\M g (T)\a) ~ T q ~2 with lower and upper bounds proved in |JSSj and [JS] respectively. Let 
us stress though that we do not claim (nor expect) that the lower bound ([3]) is optimal. 

The proof and above heuristics are based on the vanishing of the Lyapunov exponent ^ E at the 
critical energies E\ = (vr/) 2 , I G N, which have been known at least since the work of Ishii |Ish] . 
The Lyapunov exponent is defined in terms of the fundamental solution T E {x,y) G S1(2,R), 
x, y G R, of the first order linear equation in R 2 associated to (JT]), by the formula 

1 E = lim ^ log(||r^(iV,0)e||), (4) 

TV— s-oo iV 

where e is an arbitrary unit vector in R 2 . The convergence is known to be almost sure |BL] . The 
next theorem summarizes our main results on the behavior of 'j E as well as the integrated density 
of states M E in the vicinity of the critical energies. The formal definition of M E is recalled in 
Section |3j 

Theorem 2 Let E(v) 7^ 0. With the positive constants 



D E(v* - E(vf) D _ (E(v) 



16E(v)E l ' T \2Ei 

one has for e > 



and 



U El ~ £ = 1--D+E 1 * + 0(e) . (6) 
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The proof of these formulas is based on perturbation theory for products of random matrices 
around random Jordan blocks, and uses the techniques of [SS] . It is well-known that — j E +tn H 
are the boundary value of a Herglotz function w z which in terms of the Weyl-Titchmarch function 
m^j ± defined below is given by w z = E(mf, ± ). It can be read off Theorem [2] that 

E(m5± 2 ) = I - D+z? + (£>_ + iB)z + 0{z*) , 

where B is some real constant. Then Af El+£ = I + -Be + O(e^). This constant B vanishes if the 
support of the distribution of the v n is strictly positive because in this case E\ is an upper band 
edge. We do not calculate B below, but believe that it is in principle possible by the techniques 
presented below (but with considerable effort because the Fokker-Planck operators on S 1 used in 
[SS] become singular). 

Formula ([6]) shows that the density of states conserves a one-dimensional van Hove singularity 
at Ei (even though the model is random). This is due to the fact that the random potential cannot 
move any eigenvalue from below E\ to above E\ and vice versa [KN] (this follows from a basic 
Sturm-Liouville arguement, see Section [3] for details). The existence of such singularities was 
already proved in [KN], but the expansion in is new, as are the results on the Lyapunov 
exponent. For the above random Kronig-Penney model with positive potentials Gredeskul and 
Pastur [GPj have analyzed the density of states at the lower band edges (this is a Lifshitz tail 
type regime), see also [KN]. 

Acknowledgements: W. K. would like to thank Stas Molchanov and Giinter Stolz for valuable 
discussions. The authors are grateful for financial support of the DFG. 



2 Basic analytical set-up 
2.1 Definition of the operator 

Let us begin by considering Hq = —d 2 as a symmetric operator on the Sobolev space Hq(M. \ Z), 
namely ip G L 2 (R) with if}', if)" G L 2 (R) and such that ip(n) = ip'(n) = for all n G Z. The 
deficiency spaces Ker(i?Q ± i) are infinite dimensional and can be explicitly calculated. These 
deficiency spaces can be seen as direct sum of 2-dimensional subspaces attached to each point 
n G Z (that is, the local deficiency indices are (2,2) at each n). We are only interested in self- 
adjoint extensions of H that are local in the sense that they only link the left and right limits 
at one point n G Z: 

^( n ±) = h m ^(^i^); ip'{ n ±) — h m if}'{n ± e) • 
The self-adjoint extension corresponding to the 5-potentials in ([1]) is given by 

«» - Co ? ) Ct,') . 

This means that the domain J)(H W ) is given by all functions in if 2 (IR \ Z) satisfying (j7|). This 
gives a precise meaning to the operator in ([T]). If one is interested in 8'- interactions of the form 
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Y^ n w ^'ni then the matrix ^ ™ j in ((7j) is replaced by ^ . The main features described 

in Theorems [I] and [2] are also valid in this case. For sake of concreteness, we stick to the model 
(PO). Let us point out that it is also possible to consider mixed 5 and 5' potentials, but then 
Theorems Q] and [2] do not hold (because the matrices describing the boundary conditions are not 
Jordan blocks in the same basis, which is an essential element in all arguments below). 



2.2 Fundamental solutions 

Now we are interested in formal solutions ip of the Schrodinger equation H u tj} = zip at a complex 
energy z (possibly not square integrable) which satisfy the boundary conditions (J7|). As this is 
a second order differential equation, it is as usual helpful to reformulate it as a system of first 
order equation for vector valued functions 

* = (fy 6 H\R,C 2 ), *(«+) = (J *(n-) , (8) 



given by 



Conversely, solving (Q with the boundary conditions (JS} gives a solution \1/ the second component 
of which solves the Schrodinger equation H^ip = zip. These solutions define the transfer matrices, 
also called fundamental solutions, by the equation 

= T*(x,y)y(y), 

as well as the requirement to be right continuous in x and y, namely Tj(x+,y+) = Tj(x,y). 
Combined with (jSJ), this implies, in particular, 

% z (n,y) = (J V Av(n-,y), y<n. 

Furthermore, we set for x < y, 

V(x,y) = (r:(y,x))-\ 
Then the concatenation identity holds for all x,y,t £ R: 

r:(y,x) = r:( y ,t)r:(t,x), (io) 

Next let us calculate the transfer matrices on an interval in R \ Z, for example [y, x], by solving 
©: 

_ z , \ A \ A) — z\\ / cos(z^(x — y)) —z^sm(z^(x — y))\ 

% z (x,y) = exp [(x-y)( = _ L ) A y) \ } K y)) . (11) 

V V u // \ z 2 sm{Z2{x — y)) cos{Z2{x — y)) J 
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Let us note that the second expression is an even function of 22 so that the choice of square 
root is irrelevant. For sake of concreteness, we will always choose the principal value though. 
As Tj(x,y) is the exponential of a matrix with vanishing trace, it has unit determinant, namely 
Tj(x,y) G S1(2,C). For real energies Bel, one, moreover, has T u f(x,y) G S1(2,R). It is useful 
to introduce the notations 

V = 7Z(n,n-l). 

According to the above 

j-z = A v n \ ( cos(>^ -z5Bin(*5)\ ^ 
\0 1 J \z~? sin^a) cos(,2 2) J 



:-d' (J "J- 



7? = (-1)' ; ? . (13) 



Therefore, at Ei = (7r/) 2 , 



2.3 Basic spectral analysis 

Here we restrict our attention first to the non-random case of the periodic Hamiltonian H v defined 
by setting v n = v > 0. As usual, the spectrum can be read off the transfer matrices T E = 7^ 
at real energy E = k 2 e R, namely E 6 cr(H v ) if and only if |Tr(T s )| < 2. The band edges are 
therefore given by the solutions of the equation 



|Tr(T B )| = 2 cos(A;) + V - sin(Jfc) 



2 . 



Thus there exist positive constants ci(v) G M + such that 

a(^) = |J [E t - c l {v),E l ] , E l = (In) 2 

In particular, the right band edges are independent of v. Thus they are also band edges of the 
random operator if the distribution of the f n 's has positive support. Moreover, at these 
energies the transfer matrices for any periodic approximant are given by 

En 
7=m+l V j 



and therefore have a trace of modulus 2. A little bit more can be said (see |KM] and [GHK]). 

Proposition 3 Consider the Kronig-Penney model with positive independent random potentials 
distributed identically according the probability distribution p with compact support in R + . Let 
v = inf supp(p). Then the almost-sure spectrum satisfies 

a(H„) = a(H v ) . 

Proof. This follows from a standard Weyl sequence argument using operators with almost 
constant potential almost equal to v. □ 

The following variant is proved similarly. 
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Proposition 4 Consider the Kronig-Penney model with independent random potentials dis- 
tributed identically according the probability distribution p with compact support containing 0. 
Then the almost-sure spectrum satisfies 

[0,oo) C a(H u ) . 



2.4 Estimates on the transfer matrices 

Here the focus is on deriving an estimate on the transfer matrices at complex energies z = E + k, 
k G C, in terms of estimates at real energies E > 0. First of all, by analyticity one has 



V = r„ E + k r 



■n j 



for some matrix R n which is uniformly bounded in n. Using the concatenation relation ( Tl0|) one 
deduces 

n 

Tj(n,m) = T E (n,m) + k ^ TJ(n,j) Rj% E (j - l,m) . 

j=m+l 

Taking the norm in the above equation, estimating the r.h.s. and taking the supremum over 
< m < n < N leads to the following lemma [SimJ. 

Lemma 5 Set 

ci = sup \\T E (n, m)\\ , c 2 = sup||i? n ||. 

0<m<n<N n 



Then, as long as \k\ c x c 2 N < 1, 



sup \\T E+K (n, m)\\ < 



0<m<n<N 



1 — re CiC 2 N 



2.5 Green functions and Weyl theory 

At some instances below formulas connecting the transfer matrices to the Green functions will 
be used. Such a link is established by Weyl theory which is reviewed, e.g., in |KS] . For this 
purpose, one cuts H w into two half-sided operators. This cut is usually done at the orgin 0, but 
as the model ([T]) has a singular potential there, it is more convenient to cut at some other point 
a G (0, 1). Then it is a basic fact that for every complex energy z with Qm^) > there are two 
unique /* ± G L 2 (R_|-) with ± (a) = 1 and solving H^f*^ = z f% ± on R±. The m- functions are 
then defined by 

<± = (G. /« J (a) • 
In terms of the transfer matrices, one has 



It is also known that the resolvent [z — HJ) 1 is an integral operator with the following jointly 
continuous kernel 

G*(x,y) = %<y, (15) 

\fZ,+( x )( m Z,+ + m Z,-)~ 1 fZ,-(y), y<x. 
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3 Lyapunov exponent and DOS at upper band edges 



The purpose of this section is to prove Theorem [3 Let us begin with the Lyapunov exponent. 
It follows from flUJ) that 

l El = , Ex = (In) 2 . 

Then fl5]) describes how the Lyapunov exponent grows as one enters the spectrum. This will be 
achieved by a controlled perturbation theory using modified Priifer variables. This technique 
also gives access to the density of states. 



3.1 Modified Priifer variables 

The basic fact motivating the use of modified Priifer variables is that the Lyapunov exponent 
defined in (@| can be calculated by 

7 S = lim El \og(\\MTj ] (N,0)M- 1 e\\) = lim E±- log(|| (MT^M' 1 ) ••• (MT* M^eW) , 

N— ¥oo iV _/V->oo iV 

(16) 

where M is an arbitrary invertible 2x2 matrix which may, moreover, depend on E. This matrix 
M can be chosen later in such a manner that the building blocks MT^M' 1 are close to some 
adequately chosen normal form. The choice of M will be made in Section 13.31 

In order to telescop the matrix product in the Lyapunov exponent further below, let us next 
introduce a random dynamical system on the unit circle. The unit circle is identified with unit 
vectors in R 2 via 



/ COS 8\ n r„ _ x 



Then there is a natural action of invertible real 2x2 matrices on the unit circle given by 

esT{e) = w^\\ • { ] 

In particular, the map St is invertible and S^ 1 = <Sy-i. With the notation u = (_), one has 

(18) 



(u\ 


T\ 


ee) 


(u\ 


T\ 


ee) 



In connection with the calculation of the Lyapunov exponent, one now has to consider the 
random dynamical system (Markov process) on the unit circle generated by the random sequence 
MT 1 f~ £ M~ 1 of transfer matrices: 

where 9$ is some initial condition and 

The 9 n are called the M- modified Priifer phases. If one sets e n = e# n , they are explicitly given 
by 

en = WffM^AV l|eo11 = 1 • (20) 
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Coming back to the Lyapunov exponent as given in (flEl) . one now has 

1 N 

^ = N D E l0g (H^X^e^ll) . (21) 

n=l 

In particular, we now have the Lyapunov exponent given by a Birkhoff sum associated to the 
random dynamical system fl20l) (which again converges almost surely so that the expectation in 
l2~Tlmay be dropped). 

3.2 The integrated density of states 

The integrated density of states is defined by 

J\T E = lim — ^{eigenvalues of H w n < E} , 

TV— >oo N 

where H^^n is the restriction H u to [0, N] with Dirichlet boundary conditions (the Dirac potentials 
at the boundaries and N vanish). It is known that the limit defining M E exists almost surely 
and is almost surely indpendent of o>. At the critical energies Ei, one has 

M El = I. 

As pointed out in |KN] this follows from Sturm-Liouville oscillation theory because at energy 
Ei and volume iV one always (for all values of the potential) has an eigenfunction with exactly 
Nl zeros so that there are exactly Nl eigenvalues below E\. Roughly stated this means that 
there is no spreading of density of states through each E\ when the distribution of the impurities 
is changed. Furthermore, Sturm-Liouville oscillation theory allows to calculate the integrated 
density of states as a rotation number (this can be done at every energy). In the vicinity of 
the critical energies, it is, moreover, possible to use the rotation number R e associated to the 
modified Priifer phases defined by 

1 1 N ~ l 

R e = - lim - E V (S e , n {B n - X ) - 0„_i) • (22) 

n=0 

In fact, one then has 

U El ' £ = I + R e . (23) 



3.3 Choice of the basis change in the elliptic regime 

Close to the critical energy Ei = (7r/) 2 , the transfer matrix given by (fl2|) is close to a Jordan 
block. Let us begin by expanding T El ~ £ in e > 0. For sake of simplicity, let us throughout 
assume that / is even. As (Ei — e)3 = tcI — + 0(e 2 ), one finds 



T 1 





1 



0(e 2 



(24) 



S 



To lowest order, namely for e = 0, the dynamics ( TT7|) of the Jordan blocks have a unique, 
common, globally attractive fixed point 9 = which, however, is not stable. Therefore the orbit 
is close to 9 = most of the time. This feature is conserved after a small perturbation so that 
the invariant distribution of the 6* n 's is to lowest order a Dirac peak. Following |DG| ISSj . it is 
a good idea to blow up the vicinity of the fixed point in order to detect the effect of disorder 
and calculate deviations of the invariant measure from the Dirac peak. The blow up will done 
by conjugating T^ l ~ £ with 

1 



Mi 



Indeed, 



m 1 r n Ei ~ £ Mf 1 



1 



1 

£2 



v„ 



•(*S 



0{e 



(25) 



\ j 

which is close to the unit matrix. Next let us center the random variable v n . Set v = E(t> ) and 
v n = v n — E(i>). Then one can split the Mi-modified transfer matrix into a deterministic and a 



random part: 



Mi T El ~ £ M[ 



1 - S^r E2V 

I i , I - r„ 

- £2 2El 1 











0(e 



(26) 



As long as ev < AttI, the deterministic summand is up to corrections of order conjugate to a 
rotation matrix R i by an angle —rje^ where 



-T)E < 



V 



v 

2Ei 



This is what we call an elliptic regime. An equivalent condition for the elliptic regime is that 
the absolute value of the trace is less than 2. In this situation it is good to conjugate ( 1261) 
again with an adequate matrix M 2 in order to bring this rotation into the normal form given by 
cos(/3) — sin(/3) N 
sin(/3) cos(/3) 



R 



The good choice turns out to be 



Mo 



-b 
b 



Indeed, setting M £ 



M £ T 



: (M £ y 



R 



M 2 Mi one finds after some algebra 

1 1 
1 



1 + 6^ 



— e ■ 



ev r . 



2E, 



(27) 

After this basis change M £ the dynamics consists to lowest order C(e2 ) of a deterministic rotation 
and a centered random perturbation which thus has a variance of order O(e). On the one hand, 
perturbation theory of the Lyapunov exponent will be based on (1271) . but on the other it also 
allows to readily calculate the modified Priifer dynamics (fT9|) perturbatively by using ( fl8l) : 



SsM = 6 ~ V£* + (sin(20)-i; 



: £2 



+ 0{e) 



It is also possible to calculate the terms of O(e), but this is not needed here. 



(2f 
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3.4 Perturbative calculation of Birkhoff sums 

For the calculation of the Lypapunov exponent and the IDOS by (I2T]) and ff22]) respectively one 
needs to evaluate Birkhoff sums of the type 

N-l 

l N (f) = E - J2 /(*») , / e C7(§ 1 ) . 

n=0 

The following result is taken from [SSJ. 
Proposition 6 For any f G C 1 (S 1 ), one has 



Proof. Because ijv(/) = c + Ijy(f ~ c) for c = J 2,r ^ /(#), one may assume that d9f(9) = 0. 
Then / has an antiderivative F G C 2 (S 1 ). Using a Taylor expansion, 



F(0 n ) = F(S e , n (9 n ^)) = F(9 n ^) + f(9 n ^)e* (sin(20 n _ 1 ) - 1) — =£= - 77 + 0(e). 

As u„ is centered and independent of n _i, taking the expectation and summing over n shows: 

I N (F) = I N (F) - e* V I N {f) + 0(e, AT 1 ) . 
Dividing by rj finishes the proof. □ 

3.5 Calculation of the Lyapunov exponent 

Here we conclude the controlled perturbation theory for the Lyapunov exponent using the formula 
(j'2ip . In view of (I2ip . it is convenient to introduce the auxiliary random variables 

ln = log(||M^- e (MTV-J) • 

Using ( 127|) . one finds 

In = ~ - 7== cos(20 n _i)ea + i^- mn(20 n _i) e 

— ~2 

+ JL S i n (2^_ 1 ) £ + -^(l-2sm(20 n _i)-cos(40 n _i))e + O(et) . (29) 

Now w n is centered and independent of 9 n _ x so that the expectation value of the terms in the 
first line vanishes. Furthermore, by Proposition [6] 

1 N 

AT— >oo VV 

n=l 



and similarly for the Birkhoff sum of cos(40 n _i). Therefore, replacing (129 j) into (12 ip shows 

7 



TV^oo AT ^ 16 
n=l 

which is already the formula for j El ~ e in Theorem |5J 
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3.6 Calculation of the IDOS 

Next let us calculate the rotation number R El ~ e defined in (I22|) by using ( 128|) . Again using that 
v n is centered, one finds R E '~ £ = —^s^ + O(e). By (1231) this shows the formula for the IDOS in 
Theorem [2j 

3.7 Choice of the basis change in the hyperbolic regime 

In this section we consider the transfer matrices T El+£ with e > 0. Of course, in (j2l|) this only 
leads to a sign change in front of e. Then the basis change M\ gives instead of ( l26l ) 

M.T^Mr 1 = {^\ £ ^ £ "f \ +v n (^f e *) + 0(el). (30) 

The deterministic part now has trace larger than 2 and is thus conjugate to a hyperbolic matrix. 
Again using M2, M 6 and rj as above, one now finds 



(31) 

The phase dynamics in this representation is given by 

S e>n (0) = 9 + e-* V sin(2#) + — ^= (sin(20) - 1) + O(e) . 

2v2f 

It has = I as stable fixed point (and 9 = as unstable one), up to errors of order 0(e). 
Replacing this information in (1311) allows to show (by proceeding as in Section 1331) that ^ El+£ = 
rjez + 0(e) which is the missing part of Theorem |2j 

4 Transport bounds 

The aim of this section is to prove Theorem [TJ 
4.1 Set-up for the lower bound 

By the results of Section 12.51 the resolvent is an integral operator with an integral kernel given 
by f|T5|) . From this follows that also the unitary groups of have an integral kernel and we will 
denote all these integral kernels with Dirac notation. Then 

17 poo f 

(a\M g {T)\a) = - J dt e~* J dx \x\ q (a\e tHut \x) (x\e~ lH ^\a) . 

The point a is used above in the construction of the resolvent, but can be chosen freely so that 
one obtains the lower bound for all a G M \ Z. Writing out the appearing matrix elements in 
spectral representation and using the identity 

dE 1 1 



dt e-T e^e"^"* 



2tt E'-E-iT^ 1 E"-E + iT~ l 
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it follows that 



{a\M q (T)\a) 





r dE 


j dx\x\i^K 






7T 



{x\(H u - E-iT- L )- L \a) 



A lower bound is obtained by restricting the energy integral over an interval [E\ — Eq, E{\ close 
to any of the critical energies: 

f f El dE 2 

(a\MJT)\a) > dx\x\ q E / — (x\(H u - E - iT- l y l \a) 

J JEt-eo^T 

dx\x\ q E / (z|(.H w -.E , / + £-«T- 1 )- 1 |a) 

Jo 

r £ ° dp 

dx \x\ q E / \Gl{ Xl a)\\ 
Jo 

where in the last equality we set z = Ei — e + % T~ l . By (TT5T) one has for x > a 

\Gl{x,a)\ 2 = \G:(a,a)\ 2 
and similarly for x < a. Let us also note 

l^(«,a)| 2 = , * , 2 - 
Replacing thus shows, still with z = E\ — e + iT~ l , that (a\M q (T)\a) is bounded below by 



(32) 



dxx q E\ \Gl{a,a)\ 2 

TV 1 



Now it may happen that the appearing matrix elements vanish for some x, but as the transfer 
matrices vary as given by ( jTTT) between the potentials, for two points x, x' within one interval of 
R \ Z satisfying x' > x + |, one has 



J TJ(x,a)l f 



> Ci 



7^(x,a)M7: z (x,a) 

2 



T*{x,a) 



m,. 



1 



l m 5,+ | 2 + 1 
" ^ lirj(a;,a)P ' 

where Ci is a uniform lower bound (for E = E\ — e G — £o, Ej\ and T > 1) on the positive 
matrix 



.4 



l) (l) +r ^ x '^y(i) (?) v (■>■'->■) ■ 
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and where in the last inequality it was also used that the transfer matrix is in Sl(2, C) so that its 
norm coincides with the norm of its inverse. A similar bound holds for x < 0. Of course, one can 
also replace x and x' by each other at the cost of another factor. Therefore the above implies 

c r N r £o ( \m z I 2 + 1 \m z I 2 + 1 \ 

(Wla) > ^Ejf & | G> , o) |^^ w + ^ 5 L_j, (33) 

where the upper boundary iV on the space variable can be chosen at convenience later on. 
Actually, both N and Bq will be adequate functions of time T. The bound ( 1331) will be the 
starting point of the conclusion of the argument below. 

Remark Let us show how one may go over to the M £ -modified transfer matrices at this point. 
First of all, M £ = M 2 M 1 satisfies ||M £ || < C 3 and || (ikf^)- 1 1| < because ||Mi|| = 1, 

\\M^ l \\ = e"3 and C 3 = ||M 2 || = ||M 2 _1 ||. Thus (a\M q (T)\a) is bounded below by 

Ca f N f £0 9 / \m?,J 2 + l \m z u 

1 dx x q E / dee I G* (a, aY 



T Jo Jo ' wv ' n \\\M £ rj(x,a)(M £ )- 1 \\ 2 \\M e Tj(—x, a)(M £ ) _1 || 2 

Replacing (1321) for \G z u (a,a)\ 2 one can now attempt to use the deterministic lower bound 

|m^ + | 2 + l \mf„_\ 2 + l . 1 |mf,J 2 + \m z ,_\ 2 + 2 . 1 



(34) 



IK,,+I 


2 


+ 




2 


+ 2 






2 + 


ml 




2 



+ 1 — > 1 1 - 1 - > 1 , (35) 

\ m t,+ + m lj -I 2 l m w,+ + m w,-r 2 l m 5,+r + l m w,-r 2 

so that either the first or the second summand in ( |34l) is bounded below by 7. One is then 
tempted to combine this with Jensen's inequality 

HpmWf) - «p(-E(iog(iiM^(x, a) (Mr'n)) . 

Now one just calculates the expectation in the exponent exactly as the Lyapunov exponent was 
calculated. However, there is a flaw in this argument: one cannot decorrelate the two bounds 
because there is no control for which to the lower bound by 7 holds. Therefore, one needs to 
show that ||7^f(x,a)|| < e x ^ z grows as most as indicated by the Lyapunov exponent 7 2 at least 
with high probability. This is shown in the next section. An alternative would be to show that 
the two terms on the l.h.s. of f )35|) are of equal size with high probability. Another way out is 
to consider the model on the half-line with adequate boundary conditions (not Dirichlet). Then 
there is only one term and arguing with Jensen as above works. o 



4.2 Large deviations 

Here it will be shown that the M £ -modified transfer matrices appearing in remain small on 
adequate scales with high probability. For this purpose, let us work first at real energies, neglect 
the difference from x and a to the integers, and also insert an arbitrary initial unit vector eo 6 M 2 . 
Working with modified Priifer variables as in Section [3] in order to telescope the matrix product 
and then using (1251) shows 

n n 

log(||M £ r^- £ (n,m)(M £ )- 1 e ||) = £ Ti = e* ^ X j + 0((n-m)e) , (36) 

j=m+l j=m+l 
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where 



2y/2vEi 



sup \Z £ (n, m)\ > A^2" 

N<m<n<N 



These Xj are martingales, namely centered random variables depending on Vi, i < j, but not on 
the future Vi, i > j. They also depend on 9 and e. Let us begin with a standard upper bound 
on the large deviations of sums of the Xj. 

Lemma 7 Let a > 0. Set Z 6 (n, m) = X^j= m +i Xj and 
fl £ N (a, e ) = |w G Q 

Then there is a constant C5 > such that 

¥(n%(a,e )) < C 5 N 2 e~ Na . 

Proof. As the distribution of the potentials v j is compactly supported, one has — c < Xj < c for 
some c > 0. Therefore the convexity of the exponential implies 

Ej(e^) < cosh(/3c) < e^ c)2 , 

where Ej is the expectation over Vj only. For A > and /3 > 0, one thus has 

p({Z £ (n,m) > A}) = p^{ e ^ e > e ^}^ < e ~ A/3 E(e /3Ze(n ' m) ) 

= e- v E„_ 1 E„(e ft ^ E ("- 1 ' m ») < e^ 2 E n ^ z ^ n ~ 1 ^) 
< e ~ x P e ^( n - m )( c P) 2 

where the last inequality is obtained iteratively. A similar bound holds for A < and (3 < 0. 
With A = N^ +a and = N~* it follows, as n < N, 

p({|Z £ (n,m)| > A^+ Q }) < 2e^° 2 e' Na . 

Summing over m and n completes the proof. □ 
Combining the lemma with (136]) leads to the following. 

Corollary 8 For u ^ Q £ N (a, e^) one has for all —N < m < n < N and some constant Cq 

log{\\M £ T l ? l - £ {n,m){M £ )- 1 e \\) < N^ +a + C 6 Ne. 

Corollary 9 There is a set Q £ N (a) satisfying 

F(Q%(a)) < 2C 5 N 2 e~ Na , (37) 

such that for u ^ Q £ N (a), Cq as above and a constant C7 one has for all —N < m < n < N 

||7^~ e (n,m)|| < C 7 £-^ exp^ +a e^ + C 6 Ne^j . (38) 
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Proof. In order to control the norms of the transfer matrices it is sufficient to control their 
action on 2 initial vectors because for any 2x2 matrix A, 

\\A\\ = sup \\Aee\\ < V2 max \\Ae e \\ . (39) 

Hence it is sufficient to prove probabilistic bounds on ( |36i) for two initial conditions, namely one 
can set £l%(ai) = fl £ N (a, eo) Ufl £ N (a, en) which hence satisfies ( 1371) and it follows from Corollary E] 
that for u (jL Q £ N (a) and — N < m < n < N one has 

\og(\\M £ T^- £ (n,m)(M £ )- l \\) < v 7 ^ (n*+° e$ + C 6 Ne^j . 

Now the result follows from the bounds ||M e || < C3 and ||(M e ) _1 || < C%e~^ derived before 
equation [34] above. □ 

Finally, we are going to combine Corollary [9] with Lemma [5j The constant C2 in that lemma 
is equal to the r.h.s. of ( )38|) as long as w ^ Q £ N (a). Therefore, as long as 

1*1 < T^-^r = JJV^~ exp(-^ + °^ " ^10^^) > (40) 
2C2C3N N 2c 3 Cn V / 

one has for oj £ ^at(«) the bound 

||7; B; ~ £+K (n,m)|| < 2Cu£-3 exp(iV^ +a e5 + Ci A^e) . (41) 

This means that the bound fl38l) transposes, up to a factor 2, to a complex neighborhood of 
energies the size of which is given by f l40l) . 

The aim is now to use the bound (jUJ) in fl33|) and for that purpose the free parameters Eq 
and iV have to be coupled to the time T. In order for ( ]4"Ij) to be of any use for all e < e , we 
first of all choose 

£o = AT-i-2" . (42) 
Then the exponential factors in both fj4*0|) and (|4ip are of order 1 and MOD becomes 



|k| < C«iV~f" 



for some constant Cg. Now the size of the balls of radius sA = N~2~ a (in complex energy) is 
much smaller than the interval of size Sq = N~ l ~ 2a that we want to cover. However, it can be 
covered by N2 intervals of size eA. Thus setting 



j=i,...,m 



one still has a set with sufficiently small probability on the complement of which ( |41|) holds 
uniformly in an adequate set of complex energies. More precisely, we have proved the following. 
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Proposition 10 Let a > 0. There exists a set satisfying 

P(MoO) < 2C 5 N%e- Na , (43) 

such that there are constants Cs,Cg so that for all e < jV~ 1-2a and 5 < CsN~^~ a one has for 
u £ ^iv(a) 

sup \\T^ l - £+l5 {n,m)\\ < C 9 e~* . 

-N<m<n<N 

We do not expect this result to be optimal, but rather that the same bound holds for a larger 
set of complex energies given by 5 < CsN" 1 ' 01 , but could not obtain such a better estimate. 

4.3 Conclusion of the proof of the lower bound 

In ( 1331) . we now choose N and eq as 

N = (C 8 T)^ , s = N~ l - 2a = (C 8 T)-irlf . 
Due to Proposition [TU] one then has 

(a\M g (T)\a) > J dx x 9 E Xsin{<x) c ( u ) J de \G' u (a, a)\ 2 (|m^ + | 2 + |m^_| 2 + 2) 



(Cg 



,2 



where Xn N {a) c is the indicator function onto the complement of Qjy(a). Now one can use the 
deterministic bound (1351) to conclude that 



9 /"( /-JV re 

(a\M q (T)\a) > — ^ / rfxx 9 / — - (1 - P(fi*(a))) . 



Thus there are constants Ciq and Cu such that 



(a|M 9 (T)|a) > C w = ^r 1 "^ 1 )^. 

If a is redefined, this proves Theorem [IJ 
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